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Abstract 

DDD 

The theory of acoustic wave scattering by many smah bodies is developed for 
bodies with impedance boundary condition. It is shown that if one embeds many 
smaU particles in a bounded domain, filled with a known material, then one can 
create a new material with the properties very different from the properties of the 
original material. Moreover, these very different properties occur although the total 
volume of the embedded small particles is negligible compared with the volume of 
the original material. 

1 Introduction 

In [To] a theory of wave scattering by many small acoustically soft particles, embedded in a 
bounded domain D, filled with a material with known refraction coefficient, a medium, is 
developed. "Acoustically soft" means that the Dirichlet condition holds on the boundary 
of the small particles. Using the general methodology, developed in [10] , we study here the 
wave scattering on impedance particles, derive a linear algebraic system for quantities 
which yield the scattered field if the number M of the embedded particles is of order 
10, and a linear integral equation for the self-consistent (effective) field in the medium, 
consisting of the medium in which many (M — ^ oo) small particles are embedded, if 
suitable physical assumptions are made, which include the following assumptions: 

a << A, a << d, (1) 
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where a is the characteristic size of a small particle, A is the wavelength in the medium, d 
is the smallest distance between any two distinct particles. We prove that the embedded 
particles create a new material whose refraction coefficient (in the limit M oo) can be 
an arbitrary desired function, although the total volume of the embedded particles tends 
to zero as M — > oo. Thus, our theory may lead to a new technology in creating materials 
with desired properties by embedding into original material many small particles with 
the number of particles per unit volume around any point x as well as their impedances 
calculated so that the resulting new material would have a desired refraction coefficient 
n{x). The embedding of the small particles can be done using nanotechnology. 

For the theory of wave scattering by small bodies, originated by Rayleigh in 1871, we 
refer to 0, [I] and [B]. 



2 Statement of the problem and its solution. 

Let D C be a bounded domain filled with a known material with refraction coefficient 
no[x), a; e M^. The scattering problem consists of finding the solution to the equation 

Cu := [V^ + k^no{x)]u = in M^ (2) 

e*'"' / 1 \ X 

u = uo + Ao{p,a) ho - , r := |x| ^ oo, /?:=-, (3) 

r \r J r 

uo = e^'"■^ (4) 

a G S*^ is the unit vector in the direction of the incident wave, /3 is the unit vector in 
the direction of the scattered wave, /c > is fixed throughout the paper, Ao(/3, a) is the 
scattering amplitude, its dependence on k is not shown since k is fixed. 
We assume that 

no(x) = 1 in £)' := M^\D, (5) 
and write the Schrodinger equation, equivalent to ([2]): 

= [V^ + A:^ - go(a;)]M = in q^ix) := k"^ [l - no{x)) = in D' . (6) 

We assume that go is a bounded piecewise-continuous in function. Problem ([3]), 
(jl]), ([6]) has a unique solution for any square-integrable real- valued go ^ or for 

complex-valued q with Im q < 0. We sketch a proof in the Appendix. 

Denote by uq the solution to the scattering problem (l2])-(l3l), and by G the correspond- 
ing Green's function, which solves the following problem: 

CG = -6{x - y) in M^, (7) 
— ikG = \ —- \ , oo, 



d\x\ \\x 
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where S{x — y) in ([7]) is the delta function. Consider now M small particles Dm, 1 < 
m < M, embedded in D, and the scattering problem: 



dV 

U=Uo + V, -— -ikV = o{-^] , |x|^oo, (10) 







•(h) ^ 









CU = in R3\f/*i,D„, (9) 

dV 
d\x\ 
dlA 

—- = KJU onSm, 1 < m < M, (11) 
d\x\ 

where is the unit exterior normal to Sm, Sm is the boundary of Dm- We assume that 
Sm is Lipschitz uniformly with respect to m. In the impedance boundary condition (|TT1) 
the parameter h^, is a constant, possibly complex- valued, such that problem fl9l)- f|TTl) has 
a unique solution. For example, this is the case if Imhm < 0. 
Let us look for this solution of the form 

M 

U = Uo + y I G{x,t)am{t)dt, (12) 



m=l 



where am should be chosen so that the boundary condition (fTTj) is satisfied. For arbi- 
trary am G L'^{Sm) the function (fT2|) solves equation Q) and satisfies conditions (fTOj) . 
Therefore, if am are found so that conditions ( fTTl) are satisfied, then (fT2ll solves problem 
(l9l)- (fTTl) . So far we did not use the assumption that Dm are small. Let us assume ([T]), 
where ^ 

a = - max diamD^, d= min dist (Dm,Di), (13) 

and A = -j^^ is the wavelength in the medium with the refraction coefficient no{x). We 
assume that nQ{x) is practically constant on the distances of the order A. 
Let us denote 

Qm:= I am{t)dt, l<m<M, (14) 



X, t) - G{x, Xm)]crm{t)dt, (15) 



J Srr 

and rewrite f[T2l) as: 

M M 



W=Wo + X]G'(x,x„)Q„ + f^ /" [G{ 

m=l m=l 



where G -Dm can be chosen arbitrary because Dm is small. One may take Xm to be 
the gravity center of Dm- The gravity center Xm may lie outside Dm if Dm is not convex, 
but Xm belongs to the convex hull of Dm- We assume that Dm are convex, but this 
assumption is not essential. If Qm does not vanish we expect that 



\G{x,Xm)Qm\ » 



[G{X, t) - G{X, Xm)]<ym{t)dt 



Srr 



if \x — Xm\ > d » a. (16) 
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We assume that (fTBi) holds for |x — x^l >> a if Qm 7^ because 

G{x,t) -G{x,x^) = Vy G{x,y)\^^^^^^^^_^^y (t-Xm), < r < 1, (17) 
so that, for a fixed A; > 0, one has 

\G{x,t)-G{x,Xm)\ = o( ^ ) «1. (18) 

Therefore the integral in (fT6|) is 0{ka\G{x,Xm)Qm\) if Qm 7^ 0, \x — Xm\ » a, and 
ka « 1. 

In this approximation we have 

M 

U{x) =Uo{x) + ''^^G{x,Xm)Qm, \x — Xm\ > d » tt Vm, (19) 

m=l 

with an error of order 0{ka + ^). 

Since the potential qo{x) is known, one may consider the functions Uo{x) and G{x, y) 
known. Therefore the scattering problem (pl)- (fTT!) is solved if one finds Qm-, 1 < m < M. 

Let us derive the equations for finding Qm-, 1 < m < M, using the boundary condition 
([TT]) . Denote the effective field, acting on the j-th particle by Ue-, where 

^^(3.) I ^O(X) + Em^i G{X, Xm)Qm, \x - Xj\ ^ tt, ^^Q) 

\ Uo{x) + Y^^^^G{x,Xm)Qm, \x - Xm\ > d » a Vm. 
Condition (fTTj) yields: 

A _ /.^.^ ^ G(x,t)a,(t)dt = - - /^.) We(x), (21) 

where x is a point on Sj and the normal derivative on 5"^ is taken from the exterior 
to Dj domain in (1211) and below. Using formulas (A. 8) and (A. 15) from the Appendix 
one may replace the function G{x,t) in the integral in equation (1211) by the function 
go^Xjt) = ^^|^_^| with the error of order 0{ka). Using the known formula for the normal 
derivative of a single-layer potential (see, e.g., p. 5] ): 



one rewrites (|2T1) as 



dlA is] 

a^(s) = A^a^-2h,T,a,+2^^-2h,Ue{s), s e Sj, (23) 
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where 

Tja,:= I go{s,t)aj{t)dt. (24) 

JSj 

It is known (see |6, p. 96]) that 

/ ds Aja = - adt. (25) 

Integrate ([23]) over Sj, use ( IT^ and and take into account that 

f ^ds= [ mjx^V,m,{x,), V, = \D,l (26) 

where we have used the smallness of Dj to replace the integral approximately by the 
expression VjMAf.{xj) and denoted by \Dj\ the volume Vj of Dj. Thus, the integration of 
over Sj yields: 

Q =-hj [ ds [ ^j(^)^^ _ u I Ue{x)ds + VjAUe{xj). (27) 
Js, is, 47r|s-t| Js^ 



Note that 



dtaAt) [ — r ^ [ dtaiit)-^ [ dt [ — r = Q,— (28) 

where \Sj\ is the surface area of Sj, 

J -= [ [ Q.= / '^At)dt, (29) 

JsjJsj Js, 

and we have replaced the function -^i^f—jri^ practically constant at the distances of 
order a, by its mean value 

\Sj\Jsj Js,4:TT\s-t\ 

i.e., by ^^j^ | . Let us estimate the order of smallness of various terms in fl27|) . The term 
h(e{x)ds ~ Ue{xj)\Sj\ = O(a^), while the term VjAUe{xj) = O(a^), if we assume that 
Ue and AUe are bounded. Therefore, we neglect the last term on the right in (!27l) . and 
obtain 

= ^MM^^_ (30) 



1 + 



In [B] p. 27], the following approximate formula is derived for the electric capacitance Dj 
of the perfect conductor with the surface Sj-. 

C, ^ (31) 
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where the conductor is placed in a medium with the dielectric constant Eo = 1- Using 
( BT]) . one may rewrite (150]) as: 

C ~ ~ C ■ 
Qj = — lA(,{xj) := —CjlA(,{xj), Cj := — . (32) 

1 + h-^\ 1 + h-^A 

When hj oo, that is, when the impedance boundary condition becomes the Dirichlet 
condition in the limit hj — oo, then one obtains from (132|) the familiar relation Qj = 
—CjUe{xj) for the total charge Qj on the surface of the perfect conductor Dj charged to 
the potential —Ue{xj). Here Ue{xj) is defined in ( 1201) and on the distances of order a the 
field lAe{x) is practically constant. 

Substitute (l32l) into (|20|) . multiply by Cj, set x = Xj, and get 

A/ 

Qj = -QWo(x,) - G{xj,Xm)CjQm, 1 < J < M. (33) 

This is a linear algebraic system for finding the unknown Qm, 1 < m < M. The matrix 
of this system is diagonally dominant if 

i^^liY. < 1- (34) 

If condition (|34|) holds, then system (|33|) can be solved by iterations: 

^ _ J2 G{x„Xm)C, = -C,Uo{x,), (35) 

and this iterative process converges at the rate of a geometric series. Therefore, system 
(|33|) is convenient for solving the scattering problem (|9|)- (|TT1) when the number of small 
particles is not very large, M = O(IO^). If this number is very large (M ~ 10^^), then 
we study the limiting behavior of Ue{x) as M oo and derive an integral equation for 
the effective field Ue in the resulting continuous medium. 

We rewrite the second line of (l20l) using formula (132|) and obtain the following repre- 
sentation for Ue{x): 

M 

Ue{x) = Uo{x) - ^ G(X, X,n)CmUe{Xm)- (36) 
m=l 

Assume now that the limiting density of the quantities Cm exists in the following sense: 
if D C -D is an arbitrary subdomain of -D, then there exists the following limit: 

/ C{x)dx= lim y Cm- (37) 

Jd M^oo 

DraCD 
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Under this assumption one can pass to the hmit M — oo in fl36p and get 

We(x)=Wo(x)- / G{x,y)C{y)UMdy. (38) 
Jd 

Applying the operator C, defined in (Ej), to and using ([7]), one gets 

£Ue - C{x)Ue = 0. (39) 
This is a Schrodinger equation with the potential 

qix) ■.= qo{x) + C{x). (40) 
The corresponding scattering amplitude is 

A(/3, a) = AoiP, «) - ^ ^o{y, -f3)C{y)Ue{y)dy, U,{y) = Ue{y, a), (41) 

where a G S"^ is the unit vector in the direction of the incident plane wave, the following 
formula (see [5, p. 25, formula(5.1.7)]) was used: 

^ik\x\ f I \ X 

G{x,y) = ——UQ{y,-a) + o[ — ] , |a;| oo, 77=0, (42) 
An\x\ \\x\J \x\ 

and Uq^x^o) is the scattering solution, i.e., the solution to problem (I2l)-(j4l) (or, which is 
the same, to problem ([3]), (jl]), ([6])). 

To calculate the scattering amplitude by formula (14T!) one has to first solve the integral 
equation ( !38l) for Ue- The amplitude Ao{(3,a) in (|4T]) is the same as in formula ([3]). 

If the number of particles is not very large, one may use formula ( |20l) and calculate 
the scattering amplitude by the formula 

1 *^ 

A(/3,a) = Ao(/5,a) + — ^Wo(x„,-/3)g„, Qrn = Qm{a), (43) 

m=l 

where the quantities are found from the linear algebraic system fl33p . 



3 Possible applications to constructing metamateri- 
als. 

Our idea for creating a material with a desired refraction coefficient is simple: such a 
material can be obtained by embedding many small particles in the original material, 
which fills the domain D. The embedding should create the desired refraction coefficient 

n{x) := 1 — k^^q{x), where q{x) = qo{x) + C{x), (44) 
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compare with formula The function C{x) is fairly general, so the new material, that 
we have obtained from the original one, which has refraction coefficient no(x), is rather 
general. 

Let us prove that the total volume of the embedded small particles is negligible 
compared with the volume of the original material in the domain D, although the effect, 
produced by these particles on the refraction coefficient, is large. Consider a unit cube in 
D, filled with the original material. If the distance between two distinct small particles is 
not less than d, then the number of these particles in this unit cube is not greater than ^, 

3 

and the total volume of the embedded small particles in this cube is O(^). If M — > cxo 

3 

and ^ — > 0, then ^ — > 0, so that the limit of the relative volume of the embedded small 
particles per unit volume of the original material in D is zero. 

How is it possible that these particles produce large effect on the refraction coefficient? 
How is it possible that C{y) ^ 0? 

The reason is simple. Let us give a simple calculations in order to explain this reason. 
Suppose, for simplicity, that the small particles in a unit volume of D are identical. Then 
the limit (1371) exists and is not zero, provided that the limit 

hm ^ = d 

exists and C ^ 0. Indeed, O(^) is the order of the number of small particles per unit 
volume, and, by equation (|32|) . 

Cm = — = 0{a), (45) 



1 + 



a 

hm \ Sm I 



because Cm = 0{a) and the quantity Cmhm\Sm\ can be made 0(1) by choosing hm 
properly. Note that the dimension [/i^] = L~^, where L stands for length, [Cm] ~ «, 
[a] = L, [Sm] ~ a^, [Sm] = L^, so the quantity ^ '^^ | is dimensionless. 
To summarize: 

One can have limA/^oo (^) = and limM^oo ^ = C ^ 0, provided that d = 0(a^/^7), 

where 7 is a bounded coefficient whose dimension is L^/^, so that the two quantities, d 
and have the same dimension L. 

Suppose that all the embedded in D small particles have the same shape and, possibly, 
different impedances hm- Let N{x) be the density of the number of these particles per 
unit volume around point x, that is, for any subdomain D (Z D the following limit exists: 

lN{x)dx= lim ^ I- (46) 

DmCD 

One can also write 

N{x)dx = X] 1, 

DmCdx 
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where dx is a small element of volume around point x G D, such that dx still contains 
many small particles. Then 



C{x) ^ Nix)C^ = (47) 



h{x)\S\ 

Here C is the electric capacitance of a perfect conductor with the shape of a single small 
particle, \S\ is the surface area of this particle, and h{x) is the boundary impedance of 
small particles around point x. If h{x) = hi{x) + ih2{x), where hi and /i2 are arbitrary 
real- valued functions, then the function C{x) is also arbitrary. In formula (H7|) the three 
function N{x) > 0, hi{x), and h2{x) can be chosen to produce a desired function C{x), 
that is, to produce the desired refraction coefficient n{x) by formula (jH]). 

Let us give the condition on /12 and qo{x) that are sufficient for the uniqueness of the 
solution to the scattering problem flTT]) . In Appendix some sufficient conditions for 
this uniqueness have been established, namely: 

lmq{x)<0, Imh = h2<0. (48) 

If C(x) = Ci(x) + iC2{x), where Ci{x) and C2{x) are real- valued functions, then 
solving (j47|) for h{x) yields: 

CC(x) , , C C2CN(x) 
h(x) = , Imh = --- — — ^ . 49 

\S\[N{x)C -C{x)] \S\[N{x)C -Ci{x)]' + Cl{x) 

Thus, if Imh < then ImC2 < 0. Inequality Imq < holds if Imqo{x) + C2{x) < 0. 
Therefore, if one wants to create a material with the desired refraction coefficient n{x), 
i.e., with a desired q{x) = k"^ — k'^n{x), Imq{x) < 0, then one starts with an arbitrary 
no{x), ImriQ = 0, i.e., with qo{x) = k"^ — k'^no{x). Given qo{x) and q{x), one finds 
C{x) = q{x) — qo{x), and then uses formula P9|) for finding N{x) and h{x) from Ci(x) 
and C2{x). The function N{x) > in (1491) is the number of small particles per unit 
volume around point x E D. We want to prove that one can choose the functions hi{x) 
and h2{x) so that N{x) is positive, C2 < 0, Ci,C2 are given functions, and (jUj) holds. 
Let 

:= 6 > 0, H := bh{x) =Hi+ iH2. (50) 
Then the first equation fH9|) implies 

CN{x) = 9l±l^ + + 1C2. (51) 
H 

Since CN{x) > 0, and N{x) > 0, the real part of the right side of (ISTI) should be positive: 

^2 ^^2 + Ci > 0, (52) 
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and the imaginary part should vanish: 



C2H1 — C1H2 



+ C2 



0. 



(53) 



Hf + HI 



Condition (|53|) holds if 



Ci 



^Hf + Hi + H, 



(54) 



Using dSl]), write ([52]) as 



H2 ■ ^ Hl + Hi 



]>0. 



(55) 



It follows from ( l55l) that H2 and C2 are of the same sign. We may be interested in the 
materials for which the solution of the scattering problem is unique. Therefore we wish 
to satisfy conditions (148|) . The argument below shows that if Imqo = and Iraq < 0, 
then < 0. 

Indeed, if Irnqo = and Imq < 0, then H2 < 0, and fl55l) implies C2 < 0. Conversely, 
if C2 < and Jmgo = 0, then implies H2 < 0. We conclude that if H2 < 0, and 
Hi,H2 satisfy dSlD, where C2 < and Ci are given, then the number N{x) of small 
particles per unit volume around point x & D, calculated by formula (15T|) . is nonnegative 
(it can vanish around some points x), so that it has physical meaning. 

If one embeds small particles of the same shape with the density of their numbers 
N{x) and chooses the boundary impedance h{x) of these particles so that the function 
bh{x) := H{x) := Hi + iH2 satisfies the condition H2{x) < 0, and if equation (I5T1) holds, 
then the material one obtains by the embedding of these small particles into D will have 
the desired refraction coefficient n{x), and the relative total volume of these particles will 
be negligible. 

Appendix 

1. Uniqueness of the solution to the problem ([3]), (jlj), (Q. 

If Ml and U2 solve ([3]), (jl]), IQ, then v := Ui — U2 solves and staifies the radiation 
condition: 



Multiply ([6]), where v replaces u, by v, the overbar stands for complex conjugate, take 
complex conjugate of with v replacing u, multiply it by v, subtract from the first 
equation, use Green's formula, and get 




(A.l) 




(A.2) 
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Use (A.l) and let r — > oo in (A. 2). This yields 

-2ilm [ qo{x)\v\'^dx + 2ik [ A{l3) dP = 0, 

where 

likr / 1 \ 

OO. 







ay 


r := \x\ 







(A.3) 



(A.4) 



If Imqo < 0, then (A.3) implies A{l3) = 0. 

This, in turn, implies that v satisfies the following relations: 



(V^ + k^)v = in D' := \ D, 



lim 



\vrds = 0. 



(A.5) 



|x|=r 



From (A.5) one concludes v = in D' (see, e.g.. Lemma 1 in |3l p.25]). This and the 
unique continuation property for the solution of the homogeneous elliptic equation ([6]) 
imply that = in M^. The proof is complete. 

2. Uniqueness of the solution to problem (l9l)- (|TTl) . 
We argue as above and get 

-2ilm [ qo{x)\v\'^dx + 2ik [ |A(/?)p(i/? - 2? V [ Imhm\v\'^ds = 0, k > 0. 
Jd Js^ Jsm 

(A.6) 

If Imhm < 0, and Imqo{x) < 0, then (A.6) implies A{P) = 0, v \s,„= 0, and this implies 
w = in ]R^\ IJm=i ^m- This implies t; = in by the argument, given at the end of 
n.l of this Appendix. 

3. Estimates of G{x, y) as \x — y\ ^ and as \x — y\ ^ (yo. 
We start with the usual integral equation for G: 



^ik\x-y\ 

G{x,y) = g{x,y) - I g{x,t)qo{t)G{t,y)dt, g{x,y) := — 

Id 4:7r\x-y\ 



Equation (A. 7) implies 

Gix,y) = g{x,y)[l + 0{\x - y\)], as \x - y\ ^ 0, 

provided that 

f \%it}\dt 
sup / I 1 < c, 

where c > stands for various constants. We assume that (A. 9) holds. For example, 
(A. 9) holds if go is bounded. We also have 

VyG{x, y) = Vyg{x, y) - [ g{x, t)qo{t)VyG{t, y)dt. (A.IO) 

Jd 



(A.7) 

(A.8) 
(A.9) 
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This implies 

VyG{x,y) = Vyg{x,y) 
Indeed, 



1 + {\x-y\'^\n 



\x - y\ 



as \x-y\^0. (A.ll) 



Vy^(x,|/) = ikg{x,y) 1 



1 



y-x 



so 



and 



\Vyg{x,y)\^0{\x-y\-''), 



ik\x — y\J \y — X 



as — y| — > 0, 



g{x,t)qo{t)WyG{t,y)dt 



' D 

Note that 



< c 



(it 



d\x- t\\t - y\ 



< cln 



1 



F - ?/l 



= go{x,y)[l + 0(A;a)], ^o(a:,?/) := 



1 



(A.12) 
(A.13) 

as \x — y\ ^ 0. 

(A.14) 

if \x-y\< a. (A.15) 



Att\x — y\ 

Now, let |x — t/| — > oo. Assume that ^/ is in a bounded domain. Then (A. 7) implies that 

1 



\G{x,y)\ + \VyG{x,y)\^0 



\x\ 



if |a:| — > oo, \y\ < c. 



(A.16) 



References 



[1] Dassios, G. , Kleinman, R. , Low frequency scattering, Oxford Univ. Press, 
Oxford, 2000. 

[2] Landau, L. ,Lifshitz, E. , Fluid mechanics, Pergamon Press, Oxford, 1984. 

[3] Ramm, A. G. , Scattering by obstacles, D.Rcidcl, Dordrecht, 1986, pp. 1-442. 

[4] Ramm, A. G. , Equations for the self-consistent field in random medium, 
Phys.Lett. A, 312, N3-4, (2003), 256-261. 

[5] Ramm, A. G. , Inverse problems, Springer, New York, 2005. 

[6] Ramm, A. G. , Wave scattering by small bodies of arbitrary shapes, World Sci. 
Publishers, Singapore, 2005. 

[7] Ramm, A. G. , Electromagnetic wave scattering by many small particles, Phys. 
Lett. A, 360, N6, (2007), 735-741. 



12 



[8] Ramm, A. G. , Distribution of particles which produces a "smart" material, Jour. 
Stat. Phys., (to appear) 

[9] Ramm, A. G. , Many-body wave scattering by small bodies, J. Math. Phys., 48, 
2 (2007). 

[10] Ramm, A. G. , Wave scattering by small particles in a medium, Phys. Lett A, 
(2007). (to appear) 



13 



